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Abstract. The dipolar model [1] has attracted much interest because it may phenomeno¬ 
logically explain the CMB hemispherical power asymmetry found in the WMAP and Planck 
data. Since such a model explicitly breaks isotropy at large angular scales it is natural 
to wonder whether it can also explain other CMB directional anomalies. Focusing on the 
low ^ alignments and assuming ACDM, we confirm that the quadrupole/octupole and the 
dipole/quadrupole/octupole alignments are anomalous with a significance up to 99.9% C.L., 
for both WMAP and Planck data. Moreover, we show for the first time that such features 
are anomalous also in the dipolar model, roughly at the same level as in ACDM. We conclude 
that the dipolar model does not provide a better fit to the data than the ACDM. 
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1 Introduction 

Cosmic microwave background, henceforth CMB, anisotropy observations (as well as other 
astrophysical and cosmological observations) can be described with just six parameters in 
the ACDM model. To date, no extension of this model has improved in a significant way the 
fit to the available data [2, 3]. It is impressive that all the huge amount of data arising from 
cosmological observations seem to suggest that such simple model is sufficient to describe 
the large scale universe we live in. However, observed features exist that are not very well 
explained by the ACDM model. This is the case of the largest CMB angular scales where so- 
called anomalies occur. These can be grossly divided in two classes: isotropic and anisotropic 
anomalies. Examples of the former are the lack of power at large angular scale [4-7], the 
lack of correlation in the two-point correlation function [8-12] and the so-called point-parity 
anomaly [7, 13-17]. In the latter we list the hemispherical power asymmetry [7, 18-24], the 
mirror-parity anomaly [7, 25-28], the cold spot [29-32] and the low i alignments [33-41], 
The significance of these anomalies is in general of the order of 2-3 a, rarely more. 

A key point is whether these anomalies can be ascribed to residual systematic contam¬ 
ination (of astrophysical or instrumental origin), or may hint to new physics. Since we now 
know that the CMB anomalies are consistently observed in both WMAP and Planck data, 
little room is left for the possibility that they are artificially created by residual systematic 
effects. The high quality level of foreground component separation performed by Planck 
[42] appears to rule out the case for residual foreground contamination unless there are un¬ 
accounted ingredient to the foreground model, see e.g. [43] for a possible candidate. The 
simplest explanation is that of statistical flukes; such line of reasoning is supported when 
properly accounting for multiplicity of tests also known as the “look-elsewhere effect” [44], 
However, the number of these features, the fact that not all of them are related one another 
in an obvious manner and their almost exclusive occurrence at large angular scales motivate 
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the quest for a (possibly unifying) explanation even if the individual statistical significance 
is not very high 1 . 

In the current paper we focus on the low t alignments, namely the unlikely alignments 
between the quadrupole and the octupole, as well as the dipole with both of the former. In 
the light of several foreground cleaned CMB maps released by both WMAP and Planck, we 
aim at assessing the statistical significance of these features. In so doing, we test not only the 
ACDM model, but also the so called dipolar model. The latter is a phenomenological model 
which has been invoked to explain the already mentioned power hemispherical asymmetry 
[7, 9, 20-22, 45]. The dipolar model [1] consists of a particular mechanism for breaking the 
isotropy on the large-angle CMB fluctuations. The model is described by: 

(n) = (l + An-p)(^f\ (n), (1.1) 

where n is the observed direction, (A T/T) mod is the observed (and modulated) CMB tem¬ 
perature fluctuations, (A T/T) iso is the usual isotropic CMB pattern, A is the amplitude of 
the dipole modulation and p is a given direction. In [22] it is found that A = 0.07 ± 0.022, 
statistically significant at ~ 3 a and the direction p is given by (l, b) = (224°, —22°) ± 22° in 
Galactic coordinates, significant at ~ 3.3 a, see also [20] for previous results. 

The paper is organized as follows. Section 2 is the bulk of this paper. In particular in 
Section 2.1 we discuss the state of the art of the CMB anomalous alignments and describe 
the used data set. In Section 2.2 we introduce the methodology employed, based on the 
multipole vectors formalism. We set forth the estimators adopted in Section 2.3 and present 
our data analysis pipeline, employed both for real data and realistic simulations in Section 
2.4. We present our results in Section 3 while Section 4 is reserved for conclusions. 

2 CMB low i alignments 

2.1 State of the art and employed data set 

The occurrence of the anomalous alignments in the large angle CMB pattern has been noted 
since the very first appearance of the WMAP data [46]. Using a different methodology, it 
was confirmed [33] that the quadrupole and the octupole are unlikely aligned in the WMAP 
ILC 1 year data (see also [47] for a similar and independent analysis). It was later shown [8] 
that the quadrupole/octupole unlikely alignment is still present in the WMAP ILC 3 year 
map at 99.6% C.L.. Moreover in the same paper a correlation between quadrupole, octupole 
and dipole was found with a significance of 99.7%. The quadrupole/octupole alignment has 
also been studied in the Planck data [7], where similar conclusions were drawn although with 
slightly lower significance. The WMAP ILC 7 and 9 year maps are analyzed in [41] where 
it is reported that the quadruple/octupole alignment occurs with probability 0.327% and 
0.511%, respectively. In the same paper, it has been pointed out that Planck and WMAP 
data are much in better agreement after the application of the Doppler boosting correction 
[48], that is, the distortion of the CMB anisotropy pattern induced by the proper motion of 
the observer with respect to CMB rest frame. 

In this paper we analyse CMB maps from both WMAP and Planck. For WMAP we 
consider three releases of ILC (Internal Linear Combination of the multi-frequency) maps 

1 Note that such significance is largely dominated by cosmic variance in the underlying ACDM model 
assumed. 
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of the CMB sky [49], namely we use WMAP ILC 5 year [50], WMAP ILC 7 year [51] and 
WMAP ILC 9 [52], See also [53] for further details about the ILC method. While for the 
Planck satellite we use two maps of the 2013 cosmological release of data [42]: SMICA, 
Spectral Matching Independent Component Analysis, [54], that implements a parametric 
approach for foreground reduction in the harmonic domain 2 , and NILC, which employs a 
spherical needlet version of the ILC algorithm [56]. 

The alignments are visually illustrated in Fig. 1 where i = 2 and t = 3 of the Planck 
SMICA map are shown as a representative case. 


smica_ell=2,fits: TEMPERATURE 



sniica_ell=3.fits: TEMPERATURE 



Figure 1. The £ = 2 (upper panel) and £ = 3 (lower panel) contributions to the Planck SMICA map. 


2 In fact we consider an inpainted SMICA map which has been produced by replacing the masked pixels 
with a constrained Gaussian realization obtained by the method described in [55]. 
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2.2 Multipoles vectors 

It is customary to expand CMB anisotropy maps into spherical harmonics. However in the 
context of multipole alignments, it is very convenient to use an alternative and completely 
equivalent representation, namely multipole (or Maxwell) vectors expansion [33, 41, 47]. The 
fundamental idea is that the information contained in each set of (complex) ag m coefficients 
for any integer m = —£, ...,£, can be recast in £ unit (real) vectors Vi and one (real) amplitude 

aim -> A^\vi...v(.. ( 2 . 1 ) 

In fact, we note that strictly speaking the term vector is improper here because we should 
rather speak of axes or directions. This happens because the association given in Eq. (2.1) 
is defined up to a “global” sign. 

The main advantage of this formalism is that it is much easier to build quantities in¬ 
variant under rotation from multipole vectors rather than from ai m . The latter is rather 
an important point because we will make use in the following of estimators based on ro¬ 
tation invariant quantities. Unfortunately no closed analytical expression for Eq. (2.1) is 
available. Therefore, numerical routines must be used to build the vectors. Further details 
and properties can be found in [33, 41, 47]. 


2.3 Estimators 


We build eight estimators, all defined in the interval [0,1] [33, 35, 37, 41]. Of these, six are 
for the quadrupole/octupole alignment: 
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and two for the dipole/quadrupole/octupole alignment: 

DQO s = ^(\q-d\ + |oi • d\ + \o 2 ■ d\ + |o 3 • d \), 

- I [(! - \q ■ d\ ) 2 + (1 - loi • d\f + (1 - 102 • d\f + (1 - 103 • d\f] 
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( 2 . 8 ) 

(2.9) 
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DQO t = 1 


In the above equations, the symbol ' denotes the unit vector, and the area vectors q and Oj 
are defined via the following vector products: 


q = <?21 X q 2 2 , 

(2.10) 

Oi = 032 X O 33, 

(2.11) 

02 = O33 X O31, 

( 2 . 12 ) 

03 = O31 X O32, 

(2.13) 


where q 2 j (with j = 1 , 2 ) represent the two multipole vectors associated to the quadrupole 
and 031 (with i = 1,2,3) represent the three multipole vectors associated to the octupole. 
The vector d represents the dipole direction which reads ( l , b ) = (263°.99,48°.26) in Galactic 
coordinates. Note the presence of the absolute values in the definition of the estimators in 
Eqs. (2.2)-(2.9) which is due to the fact that multipole vectors define directions, i.e. they are 
headless vectors, see Section 2.2. 

The estimators introduced in Eqs. (2.2)-(2.9) can be divided in “S” and “T” statistics as 
denoted by the labels. They measure “distance” from a situation of complete misalignment, 
i.e. orthogonality, which is associated to zero in both cases, whereas complete alignment, i.e. 
parallelism, is represented by the value 1. However, the “S” estimators weight the cosine 
contributions from the scalar product linearly while the “T” estimators weight it quadrat- 
ically. Note that in principle these two sets do contain different statistical information but 
we anticipate that they provide very similar results [40]. 

2.4 Simulations pipeline and observed data analysis 

We perform 10 5 Monte Carlo simulations, extracting a^ m coefficients from the Planck 2013 
ACDM fiducial model 3 . For each realization, we transform to multipole vectors employing 
the publicly available code written by Copi et al. [33], whose use is acknowledged here 4 . 
Then, for each of the performed realizations we compute the eight estimators defined in 
Eqs. (2.2)-(2.9). We therefore can build the empirical distributions of the estimators in the 
ACDM model, see green histograms in Fig. 2. For the dipolar model, our pipeline flows in 
a similar way. The only difference is that once the a^ m are drawn, we transform them to a 
real space map, i.e. AT/T\i so , and use Eq. (1.1) to compute AT/T\ mod . We then go back to 
harmonic space, i.e. 

AT/T\ mod ^aZ d , 

and use these to compute the multipole vectors. Once this is repeated 10 5 times, we 
can build the eight empirical distributions of the considered estimators in the dipolar model, 
see the red histograms in Fig. 2. 

Of course the same estimators are evaluated for five observed CMB maps, see Section 
2.1. These values are represented by the vertical lines in Fig. 2: WMAP ILC 5 in blue, WMAP 
ILC 7 in pink, WMAP ILC 9 in balck, Planck 2013 NILC in cyan and Planck 2013 SMICA 
in magenta. In fact before evaluating these numbers, we have applied a “boost correction” to 
the observed a^ m coefficients. This is necessary because the observed quadrupole is slightly 
affected by the motion of the satellite with respect to the CMB rest frame. The details of 
this correction for every multipole i are given in the next subsection. 

3 We have tested that the particular model chosen is irrelevant. 

4 See http://www.phys.cwru.edu/projects/mpvectors/ 
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2.4.1 Boost correction 


It is possible to show, see e.g. [57, 58], that the spherical harmonic coefficients, afp. observed 
in the CMB rest frame (hereafter S cm b ) are related to the spherical harmonic coefficients, 
a( m defined in a frame S which is moving in the z direction at velocity v with respect to 
S cm b, hr the following way 

OO 

a 'lm = ^ a Vrn.I™e( V ) ) ( 2 - 14 ) 

£>=0 

where no sum on m is understood and where the ipJv) is defined as 

I»] = [ +1 27T y /±E^pp[x)if dx, (2.15) 

J -1 1 + vx \l + vxj 

with the Pp functions defined through the Legendre polynomial P"' as 


P 


2l + l(l-m)! 

4t r (£+1)! * 


(2.16) 


In fact we need to invert Eq. (2.14) and “deboost” the WMAP and Planck observations. 
This can be done using the following orthonormality relation 



£' 

(2.17) 

and considering that 

mv) = m- v ). 

(2.18) 

Therefore one finds 

a fm = a e'ml™e(— v ) ■ 

(2.19) 


l' 


In practice, only 1 = 2 has to be corrected by this kinematic term. For this multipole, the 
typical correction is roughly around 10 — 30%. For i > 3 this effect is completely negligible. 
For the octupole the maximum deviation is computed to be of the order of 0.1%. See 
Appendix A where explicit values are reported. 


3 Results 

Our results are shown in Fig. 2. We evaluate the level of anomaly comparing the histograms 
with the observed values, i.e. the vertical bars in Fig. 2. We consider both the ACDM and 
dipolar, and for each analyzed CMB map, i.e. WMAP ILC 5, WMAP ILC 7, WMAP ILC 
9, Planck 2013 NILC and Planck 2013 SMICA. 

At the price of a slight inaccuracy in terminology, we define the probability to exceed, 
henceforth PTE, as the number of the simulated counts that have the value of the considered 
estimator smaller that the observed value. These values are reported in Table 1 and the 
PTEs are provided in Table 2 and in Table 3. 
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Figure 2. S statistic for the upper row and T statistic for the lower row. Green histograms for the 
empirical distribution of the considered estimators in ACDM and red for the dipolar model. From left 
to right we consider S, S 23, S23 and DQOs in the first row and similarly T, T23, T23 and DQOt in 
the second row. Vertical lines are for the observed estimators (already boost-corrected): WMAP ILC 
5 in blue, WMAP ILC 7 in pink, WMAP ILC 9 in balck, Planck 2013 NILC in cyan and Planck 2013 
SMICA in magenta. In each panel we show the counts in the y-axis and the estimator in the x-axis. 


Table 1. Values of the estimators extracted from the WMAP and Planck CMB maps. 


Estimator 

WMAP ILC 5 yr 

WMAP ILC 7 yr 

WMAP ILC 9 yr 

Planck SMICA 

Planck NILC 

S 

0.799 

0.804 

0.807 

0.794 

0.804 

T 

0.959 

0.962 

0.963 

0.956 

0.962 

S23 

0.776 

0.783 

0.788 

0.718 

0.697 

T23 

0.949 

0.953 

0.955 

0.919 

0.908 

§23 

0.869 

0.877 

0.884 

0.859 

0.877 

T23 

0.982 

0.984 

0.986 

0.979 

0.985 

DQOs 

0.789 

0.792 

0.799 

0.774 

0.776 

DQO-p 

0.940 

0.943 

0.946 

0.936 

0.944 


A few comments are in order. First, the empirical histograms for ACDM and dipolar 
model are very similar. This means it is not easy to distinguish between the two models 
on basis of the observed alignments. Second, all vertical lines are very close to each other. 
This means that at large angular scale in temperature the CMB maps obtained with two 
different experiments and with three different methods are very similar in terms of phases. 
Third, all vertical bars, for all the considered estimators, stand in the right hand part of the 
histograms. This means that data tend to show alignments of the considered low multipoles. 
The significance of these alignments is in general larger than 99%, with few cases at the level 
of 98 — 99%, and can be as large as 99.9% in selected cases. 
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Table 2. Percentage of anomaly for the quadrupole/octupole alignment, for all analysed estimators 
(5, T, 523, T23, 523 and f 23) for the WMAP data (WMAP ILC 5 yr, WMAP ILC 7 yr and WMAP 
ILC 9 yr) and for the Planck data (Planck SMICA and Planck NILC). 


Estimator 

WMAP 

ACDM 

ILC 5 yr 

Dipolar 

WMAP 

ACDM 

ILC 7 yr 

Dipolar 

WMAP 

ACDM 

ILC 9 yr 

Dipolar 

Planck 

ACDM 

SMICA 

Dipolar 

Planck 

ACDM 

NILC 

Dipolar 

S 

99.647 

99.640 

99.701 

99.701 

99.750 

99.731 

99.581 

99.578 

99.704 

99.707 

T 

99.828 

99.832 

99.856 

99.866 

99.873 

99.880 

99.775 

99.769 

99.856 

99.866 

S23 

99.722 

99.724 

99.793 

99.791 

99.838 

99.830 

98.649 

98.606 

97.990 

97.951 

T23 

99.863 

99.868 

99.892 

99.891 

99.905 

99.906 

99.217 

99.207 

98.861 

98.833 

§23 

98.355 

98.308 

98.569 

98.539 

98.689 

98.676 

98.128 

98.089 

98.550 

98.523 

T23 

98.654 

98.646 

98.839 

98.802 

98.901 

98.881 

98.420 

98.379 

98.839 

98.806 


Table 3. Percentage of anomaly for the dipole/quadrupole/octupole alignment, for all analysed 
estimators (DQO_S and DQO_T) for the WMAP data (WMAP ILC 5 yr, WMAP ILC 7 yr and 
WMAP ILC 9 yr) and for the Planck data (Planck SMICA and Planck NILC). 



WMAP 

ILC 5 yr 

WMAP 

ILC 7 yr 

WMAP 

ILC 9 yr 

Planck 

SMICA 

Planck 

NILC 

Estimator 

ACDM 

Dipolar 

ACDM 

Dipolar 

ACDM 

Dipolar 

ACDM 

Dipolar 

ACDM 

Dipolar 

DQO.S 

99.803 

99.796 

99.829 

99.823 

99.872 

99.865 

99.672 

99.662 

99.687 

99.681 

DQO_T 

99.776 

99.779 

99.810 

99.808 

99.859 

99.851 

99.725 

99.728 

99.825 

99.823 


4 Conclusion 

We have tested the CMB quadrupole/octupole and dipole/quadrupole/octupole alignments 
for several foreground cleaned products for both WMAP (5, 7 and 9 year data) as well as 
Planck 2013 data. Specifically, we have considered the WMAP ILC products for the several 
year releases and Planck NILC and SMICA maps. We have used a total of eight estimators 
based on the multipole vector formalism, two for the dipole/quadrupole/octupole and six for 
quadrupole/octupole alignments. All these estimators are supported by a large Monte Carlo 
of 10 5 independent maps. We report that all the data combinations and all the estimators 
we have tested exhibit anomalous alignments for both combinations of multipoles considered, 
typically at the 98%-99% level, and up to 99.9% in selected cases. The consistent pattern 
for the alignments observed in both WMAP and Planck strongly disfavours an origin of 
the effect related to unaccounted instrumental systematics. The wide frequency leverage 
of the Planck data (30 to 353 GHz), weakens considerably the case for residual foreground 
emission. The fact that we find consistent results also among different foreground separation 
procedures (SMICA, NILC and WMAP’s ILC) makes this conclusion stronger. We have also 
investigated the possibility that the phenomenological dipolar model may provide a better 
framework for the existence of the observed alignments with respect to plain ACDM. This 
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possibility is, in principle, intriguing because the dipolar model has gathered some success 
in accounting for other anomalies, e.g. the hemispherical asymmetry. We report negative 
findings on this last issue: the dipolar model does not seem to be able to accomodate for the 
existence of anomalies significantly better than ACDM. 
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A Impact of the boost correction 

In Table 4 and 5 we report the a^ m for quadrupole and octupole without and with de-boosting 
correction respectively. 


Table 4. agm for £ = 2 and i = 3 (no correction applied). Units: /iK. 




a 20 

a 21 

a 22 

a 30 

a 31 

a 32 

a 33 

WMAP ILC 

5 yr 

12.350 

-1.087+6.069i 

-14.211-17.858i 

-6.449 

-12.733+2.443i 

22.019+0.698i 

-11.813+33.393i 

WMAP ILC 

7 yr 

11.771 

-0.771+6.215i 

-14.120-17.941i 

-6.479 

-12.191+2.026i 

21.999+0.591i 

-11.709+33.554i 

WMAP ILC 

9 yr 

12.563 

-1.727+6.209i 

-13.846-18.017i 

-6.844 

-11.271+1.581i 

21.857+0.535i 

-12.060+32.853i 

Planck 2013 SMICA 

13.089 

-1.530+2.497i 

-15.503-17.091i 

-5.959 

-12.841+1.671i 

22.086+1.670i 

-12.465+29.402i 

Planck 2013 ! 

NILC 

13.512 

-1.375+1.722i 

-13.564-16.325i 

-6.117 

-9.547+1.896i 

22.242+1.875i 

-12.914+28.340i 
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